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ABSTRACT 


In this report we present a time-varying sliding mode control (TV-SMC) technique for 
reusable launch vehicle (RLV) attitude control in ascent and entry flight phases. In ascent flight the 
guidance commands Euler roll, pitch and yaw angles, and in entry flight it commands the 
aerodynamic angles of bank, attack and sideslip. The controller employs a body rate inner loop and 
the attitude outer loop, which are separated in time-scale by the singular perturbation principle. The 
novelty of the TVSMC is that both the sliding surface and the boundary layer dynamics can be 
varied in real time using the PD-eigenvalue assignment technique. This salient feature is used to 
cope with control command saturation and integrator windup in the presence of severe disturbance 
or control effector failure, which enhances the robustness and fault tolerance of the controller. The 
TV-SMC is developed and tuned up for the X-33 sab-orbital technology demonstration vehicle in 
launch and re-entry modes. A variety of nominal, dispersion and failure scenarios have tested via 
high fidelity 6DOF simulations using MAVERIC/SL IM simulation software. 


1. INTRODl CTION 

Flight control of both current and future reusable launch vehicles (RLV) in ascent and descent 
modes involves attitude maneuvering through a wide range of flight conditions, wind disturbances, 
and plant uncertainties including aerodynamic surfac es and engine failures. The baseline RLV flight 
control system that was designed for the X-33 technology demonstration sub-orbital launch vehicle 
employs a variable structure PID control law 1 with gain scheduling. This requires four gains per 
channel that are looked up from a table as a function of relative velocity. Depending on the flight 
trajectory, each gain table can have as many as 25 values, so potentially 300 gain values must be 
stored in the on board computer for nominal flight. In case of an engine failure, or Power Pack Out 
(PPO) alternate sets of gain tables are used, depending on the flight time when the failure occurred. 
Provisions are made for 25 possible PPO times, or 25 sets of PID tables. This amounts to another 
7500 values to be stored, or a total of 7800 values to provide gains for the nominal and engine failed 
cases. The reason for so many gain tables is because the control system design relies on linear 
analysis and perturbation theory at specific design points along the trajectory. This method is well 
established and has been used in many launch vehicle control system designs. Robust control is 
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PROJECT OBJECTIVES 


The proposed research is dedicated to improving sliding mode controller for reusable launch 
vehicles (RLV) of second generation to address aerodynamic surface failures and thrust 
deficiencies. The proposed research effort is an extension of the research done under the project 
titled “Advanced Guidance and Control for the RLV of Second Generation.” This proposal 
contributes to risk reduction area associated with reduced controllability of RLV in presence of 
control effector failure and RLV model uncertainties. Availability of this type of control system 
would improve safety of the RLV, reduce operation cost for the RLV by reducing the pre-flight 
effort for each flight significantly. Increasing reliability/life of the entire system will be also 
achieved. 

Flight control of the RLV of second generation involves large attitude maneuvers through a 
wide range of flight conditions from lift-off to entiy. Improving recovery of tracking performance 
of the RLV from the aerodynamic surface failures and thrust deficiencies is achieved via re- 
configurable/adaptive time varying multiple loop continuous sliding mode controller (TV-SMC) 
design with a direct adaptation taking into account the torque command limitations. Providing a 
tracking performance recovery, the re-configurable adaptive time-varying continuous sliding mode 
controller will achieve robust, high-accuracy tracking of guidance trajectories for large attitude 
maneuvers through a wide range of flight conditions in presence of aerodynamic surface failures 
and thrust deficiencies and model uncertainties. The designed controller also will automatically 
adjust to changing specifications, such as mass of payload and target orbit, and the operating 
environment, such as atmospheric perturbations and interconnection perturbations from the other 
subsystems of the RLV. The TV-SMC will be developed and tuned up for the X-33 sub-orbital 
technology demonstration vehicle in launch and re-entry modes via high fidelity 6DOF simulations 
using MAVERIC/SLIM simulation software. A variety of failure scenarios will be tested. 
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ensured as long as the vehicle performance and operating conditions are relatively close to the 
design points. A robust flight control algorithm that would accommodate different trajectories and 
aerodynamic surface and engine failures without gain scheduling would be an improvement over 
the RLV current flight control technology. Sliding Mode Controller (SMC) is an attractive robust 
control algorithm for the RLV ascent and descent flight controller designs because of its inherent 
insensitivity and robustness to plant uncertainties and external disturbances 2 ' 4 . Such a robust 
controller 5 ' 6 would reduce risk and drastically decrease the amount of time spent in pre-flight 
analysis, thus reducing cost. 

In this work we present a time-varying sliding mode controller (TV-SMC) design technique. 
The RLV fixed-gain (FG) SMC two-loop structure, which is developed in the works 5 ' 6 , is employed 
in this work for the TV-SMC design. 

In the outer loop, the kinematics equation of angular motion is used with the outer loop 
SMC to generate the angular rate profiles as virtual control inputs to the inner loop. In the inner 
loop, a suitable inner loop SMC is designed so that the commanded angular rate profiles are 
tracked. The inner loop SMC produces roll, pitch and yaw torque commands, which are allocated 
into end-effector deflection commands. Multiple (ime scaling (multiple-scale) is defined as the 
time-constant separation between the two loops. That is, the inner loop compensated dynamics is 
designed to be faster than the outer loop dynamics. The resulting multiple-scale two-loop SMC, 
with optimal torque allocation, causes the angular rate and the Euler angle tracking errors to be 
constrained to linear de-coupled homogeneous vector valued differential equations with desired 
eigenvalues. 

In this paper the TV-SMC design technique is presented for the attitude controller (autopilot) for 
RLV. Any partial or complete failure of control actuators and effectors will be inferred from 
saturation of one or more commanded control signals generated by the controller. The saturation 
causes a reduction in the effective gain, or bandwidth of the feedback loop. The controller cannot 
tell nor does it care whether control command saturation is due to wind gust or a stuck effector, 
because it only computes the required torque to eliminate the tracking error. A truly adaptive 
controller should be able to do this even if the vehicle has lost significant control authority due to 
control effector failure. It is up to the combination of the control allocation and the control effectors 
to realize the required torque. In order to maintain stability, the bandwidth of the nominal (reduced- 
order) system will be reduced accordingly using a time-varying bandwidth PD-eigenstructure 
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assignment technique 7,8 . The presented TV-SMC fault-tolerant technique automatically handles 
momentary saturations and integrator windup caused by excessive disturbances, guidance command 
or dispersions under normal vehicle conditions. 

2. RLV ATTITUDE CONTROL PROBLEM 

The rigid body equations of motion for an RLV is given by ( 1 )-(5) 

(J 0 +AJ)<y = -Q(J 0 +AJ)ty+(I + E)T + d 

y = R (y)co 

0 -co i co 2 

Q = o) } 0 -<y, (3) 

- co l cu, 0 

where the rotational matrix R(y) for ascent is given by 

1 tan 6 sin <p tan 6 cos (p 

R(y) = 0 cos <p -sin (p y = [<P 9 vY > (4) 

sin <p cos (p 

cos 6 cos Q 

and the R{y) for entry 

cos or 0 sin or 

R(y)= 0 10 , y = [o h or j3] T . (5) 

sin or 0 -cos or 

The control problem for the RLV in ascent and descent modes is to determine the control torque 

command vector T such that the commanded orientation angle profiles y c are robustly 

asymptotically followed in the presence of bounded disturbance torque d, the RLV inertia 
variations AJ and aerodynamic surface and engine failures that are described by the uncertain 

matrix AD(.) , E = AD(.)D\ D*(.) = D(.) t [d(.)D(.) t ]"', D( ,)is a nominal sensitivity matrix. 

3. SUMMARY OF THE SMOOTH MULTIPLE-LOOP FIXED-GAIN SLIDING MODE 

CONTROL 

The smooth multiple-loop fixed-gain sliding mode controller is designed as follows 5,6 : 


( 1 ) 

( 2 ) 
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( 6 ) 


The outer loop smooth fixed-gain SMC generates body rate commands 
CQ c =R _1 (x)[f c +K,?'J + R“ 1 (»K 2 <7 

t 

G = Y e +K i jr e dr, a<= R 3 , y e =y c -y 


( 7 ) 


or 


<o c =R"'(y) 


l 

y c +(K l+ K 2 )y e +K 2 K l jy € dr 


( 8 ) 


that provide the following outer loop compensated dynamics (given the (O c profile is tracked 
perfectly in the inner loop): 

f t +(K l +K 2 )f t +(K 2 K l )y t =0. (9) 

• The inner loop smooth fixed-gain SMC generates control torque commands 

T = J o&c + Jo K 3&e +&J 0 co+J 0 K 4 s (10) 


f 

s = + K 3 jco e dr , seR 3 , co e =co -o) 


( 11 ) 


or 


T = J 0 ^c +^J 0 « C +(J 0 K 3 -QJ 0 -i- J 0 K 4 )ty e + J 0 K 4 K 3 jo) e dr ( 12 ) 

0 

that provide the following inner loop compensated dynamics (given the T profile is allocated 
perfectly into commands to actuator deflections): 

cd e +(K 3 +K 4 )y e +(K 4 K,)y e =0 (13) 

• Writing inner and outer loop tracking error equations in a damping factor/natural 
frequency format 

x + 2^O) n x + Q) 2 „x = 0 ( 14 ) 

it is easy to calculate elements of diagonal mairices K, K 2 ,K 3 and K 4 to provide given 

damping factors (usually all are equal to 1.1 or so) and natural frequencies that provide for a 
sufficient time-scale separation between the control loops to eqs. (9) and (13) for inner and outer 
loop compensated tracking dynamics. 

• Commands to the actuator deflections can be calculated as follows 


6 



( 15 ) 


T=J 0 <u, .+ru 0 w,+(J fl K,-m 0 +p 2 J 0 Kj< > e +p 2 J 0 K 4 K,ja>,</r 

0 

where S c0 is a bias command that should trim the RLV, D(.) is a sensitivity matrix (for a 

feedforward control allocation case). 

4. SMOOTH MULTIPLE-LOOP TIME-VARYING SLIDING MODE CONTROL 

Motivations. The elements of the command torque vector (12) are limited by physical abilities 
of the RLV flight control system that implies the following limits to elements of the vector T : 

a t QZT t Zb, (.), < = U. (16) 

a^^T^bX), i = U 07) 

at 

where a, (.), b t (.), a] (.), b] (.) depend on the RLV current flying conditions including Mach number, 
dynamic pressure and trim conditions. Designing the time-varying SMC inequalities (17) could be 
originally out of consideration. 

The components of the vector T that depend on body rate tracking profile, w c , and its 
derivative, w c , usually have the largest amplitudes among the others. So, the following idea is 
proposed: if the inequalities (16) and/or (17) are about to be violated we have to start reducing 
magnitudes of the terms co c and cb c to prevent actuators from saturation. It will be achieved via real 
time adjusting bandwidths of inner and outer flight control loops. This procedure requires 
application of time-varying linear control technique 7,8 that is incorporated into the time-varying 
sliding mode controller design. 

The following outer, co c , and inner, T, continuous, time-varying sliding mode control (TV- 
SMC) laws are designed 

0) c =R-\r)[7c +K,(0r e ] + R" l (T)K 2 (0o- (18) 

T = J 0 ty c + J 0 K 3 (0<y e +QJ 0 O) + J 0 K A (t)s (19) 

where the outer, a , and inner, 5 , sliding quantities are identified 

t 

a = r e +(K l (t)-K 2 (ty l k i (t))jr^T,cx€R i , y e =y c -y (20) 

0 
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( 21 ) 


S = 0) e +(k,(/)-K 4 (0 l K 3 (0)j6> e </T > J€ R 3 , (O e =Q) C -CO 

0 

Remark. One can observe an apparent difference between the eqs. (20), (21) and eqs. (7), (11). 
These equations become identical if K, = const and K 3 - const 

The diagonal elements of the matrices K 2 (t) at d K 4 (0 determines the time-varying thickness 
of the boundary layers of the outer and inner loops, respectively. Once K 2 (t) and K 4 (0 are fixed, 
K,(0 and K 3 (0 define the dynamics of the sliding modes of the outer and inner loops, 

respectively. When these controller gains are chosen in a certain manner as defined below, they can 
be used to define the desired time-varying bandwidths of the outer and inner loops. These 
bandwidths can be adjusted in real time to prevent actuators from saturation. The corresponding 
bandwidths or natural frequency of the compensated loop dynamics are to be decreased if the 
actuator command is close to saturation and are to be increased otherwise. 

5. CONTROLLER DESIGN STABILITY ASSESSMENT: NOMINAL CASE 

In the nominal case where E = 0, AJ = 0, d ~ 0, the mathematical model of the RLV in eqs. 
(l)-(5) is completely known. For the outer loop, substituting eqs. (18) and (20) into eq. (2) and 
differentiating eq. (2) we obtained the compensated outer loop error dynamics 

Ye + (K , (0 + K 2 (t))y e + (K , (OK 2 (0 - K , (Ofc = Rffl, , y e = J y e (v)dr (22) 

0 

Assume that the transient response in the inner, "faster", loop is stable and co = co c in eqs. (2) and 

(18). Then the tracking error a> e is stabilized in the inner loop at zero. Assuming ft5 c =0 the 
equation (22) can be rewritten in a homogeneous format 

y e + (K , (0 + K 2 (0)n + (K , (OK 2 (0 - K , (t))y e = 0 (23) 

Since the matrices K,ft) and K 2 (t) are diagonal, eq. (23) can be rewritten in a scalar format 

h + (*1/ (0 + *2, (0 )Yei + (*„ (0*2/ (0 - K (ofe,- =0, V/ = U (24) 

The series D-eigenvalues of LTV differential equations (24) can be identified as follows 

K (0 = (0, K , (0 = * 2 i (0 , Vj = U . (25) 

In order to assign values for the bandwidth (cut-off frequency or natural frequency) and damping 
factor for the compensated dynamics of the outer loop eqs. (24) are rewritten in a format 
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hi +a>l(t)Y tl =0, £ >1 Vi' = 1,3 . 


(26) 


Yei 


( 


2^« n ,(0 


<(0 " 

"m(0> 


Then the so-called "PD-eigenvalue" can be identified as follows: 


A,(0=(-6 

A 2 (0=(-^ -V^. 2 “ 1 ) ft, nl(0 


(27) 

(28) 


Note that for time-varying (O ni (/) , the relations between the sign of the radicals and the index of A, 
must be maintained. Then the SD-eigenvalues are related to the PD-eigenvalues as follows 

A (0 = A,<0 - (-6 (') = -MO .Vi = U (29) 

or 

*„(0 = ({,- i /* J ->)»-(0 < 30 > 


and 


MO = A*(0 + ^ - f- f, - ^ ' K <0 + ^ = -MO (31) 

a ni (0 v 7 co m (t) 

or 

*2,(0 = ff,+7^V,(0-^ (32) 

V J ^ni(0 

Thus, given and (O m ( t ) the corresponding time-varying coefficients in the sliding surface (8) can 
be computed using eqs. (30) and (32). 


Remarks. 

1. The outer loop dynamics are obviously decoupled and the stability is guranteed by the PD- 
eigenvalues for £, >1 and co ni {t) > 0, V/ > : 0 , which can be achieved by the choice of 
diagonal time- varying matrices K, (t), K 2 (t) according to (30) and (32). 

2. The gain matrices K,(t), K,(r) must be selected in accrdance with (30) and (32) with £ > 1 
in order to to make <7- dynamics faster then y e -dynamics. The time scale separation increases 
with . Note that (30) and (32) are for constant £. only, and as such £ should not be 
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adjusted in real time. The PD-eigenvalue synthesis formulas for time-varying £.(/) is 
available, but it is much more complex. 

The inner loop compensated error dynamics are described by the eqs. (1), (3), (19) and (21), 
which can be combined as 


co e + (K 3 (0 + K 4 (i t))cb e + (aT 3 (t)K 4 (0 - a: 3 (r)K = 0 , co e = 'jco e (T)dT 

0 

Eq. (33) can be rewritten in a scalar "damping factor-natural frequency" format V/ = f3 


(33) 


(O oi + 


2|A, (0 


^ni(0 

®i(0; 


^+^(0^=0, £>i. 


(34) 


Eqs. (33) and (34) is a full anology to eqs. (24) and (26). So, the following coefficients can be 
identified by analogy 


^ii (0 ~ {jfi VS (0 (35) 

( 36) 

Now, given and <y ni (0 the corresponding time-varying coefficients in the sliding surface (21) 
can be computed concerning eqs. (18) and (19) Vi = 1,3 . 


Remark. The same comments for K,(/),K 2 (0 following (30) and (32) are applicable to 
K 3(0,K 4 (0 in relation to (35) and (36), and the s- djmamics and co e -dynamics. 


6. STABILIT Y ASSESSMENT: PERTURBED CASF 

In the perturbed case assume E * 0, AJ * 0, d*0, unknown but bounded, i.e. the math 
model of the RLV in eqs. (l)-(5) is not completely known. For the outer loop, again assume that the 
transient in the inner "faster" loop is stable and a) = co in eqs. (2) and (18). Then the perturbed case 

coinsides with the nominal case, since eqs. (2), (4) and (5) do not contain uncertaintes and 
disturbances. 

The inner loop error compensated dynamics are described by the eqs. (1), (3), (19) and (21). In 
the perturbed case this is 
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( 37 ) 


& e +(K } (0 + k 4 (t))S) e + (a: 3 (t)K 4 (t I - a: 3 (ofe = F(B e , a > e , 4 , t) 
where F{(o e ,(b c ,(d e ,t) depends on perturbation teims. The closed-loop stability of the perturbed 

system can be justified by the well-known results from the Lyapunov analysis 9 . 

Theorem. Given a perturbed linear time-varying dynamic system 

x = A(t)x + F(x,t) (38) 

where ||/ r (x,t)|| ^ 8 is a bounded, nonvanishing perturbation, i.e. F(0,t ) ^ 0, and the matrix A(t) 
is Hurwitz, i.e. there exists a pair of positive definite matrices (P(t), Q(tj) satisfying the Lyapunov 


equation 


P(t) + P(t)A(t) + A T (t)P(t) = -Q(t ) 

(39) 

such that 


c \ ||*(0|| 2 ^ X T (t)P(t)x(t) < c 2 |x(0| 2 

(40) 


•* T (OQ(t)x(t) < -c 3 ||x(0|| 2 

(41) 


t r 

VI 

g 

H 

* 

(42) 

for all t > t 0 , 

where ||v|| = Vv 7 v . For any 0 < Q < 1 , define 



\(t»= 2c \Fe 

c 4 y c 2 

(43) 


If ||T(x,0|| <8<8 0 , then every trajectory of the perturbed system is exponentially, ultimately 
bounded in the sense that for any r> 0, and ||x(t 0 )|| < r Ik, there exists a t x >t 0 such that 

||x(/,/o)||^%(to)|| e '' l( '" ,o) , t 0 <t<t ] , (44) 

||x(r,r 0 )||<^ t>t x , (45) 

where 


k = 



(1 ~g)c 3 

2c 2 


, b- 


8 

8 0 (0) r 


7. SIMULATION 


( 46 ) 


The results of following simulation tests performe d on Slim 1.2 are presented: 
1 Nominal flight evaluation. 
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2 Dispersion case evaluation. 

3 Actuator-failed flight evaluation 

Testl. Nominal Flight Evaluation 

The goals and implementation conditions of the nominal flight evaluation test are 

• Capability for generating control commands that assure RLV stability and high tracking 
accuracy is to be demonstrated. 

• The aerodynamic model and actuator model is based on X-33 vehicle data. 

Roll pitch and yaw angle tracking in ascent shown in Figures 1, 2 and 3 demonstrate high tracking 
accuracy. 

The aerodynamic surface deflections that are shown in Figures 4 and 5 do not hit the limits. 

Body torque ratios that are shown in Figure 6 demonstrate the fact that commanded torques are far 
from saturations, that’s why the gain adaptation algorithm is not activated. 

Angle of attack, bank angle and sideslip angle tracking in descent shown in Figures 7 - 9 
demonstrate high tracking accuracy. 

The aerodynamic surface deflections that are shown in Figures 10 and 1 1 do not hit the limits. The 
RCS propellant used in descent is shown in Figure 1 2. 

The corresponding body torque ratios are far from limits, and the gain adaptation algorithm is not 
activated in descent nominal mode. 
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Fig. 9 Sideslip angle tracking (degrees) 




Fig. 10 Elevon deflections (degree) 
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Summary’ of the Test l Results 

• The nominal mission trajectory was successfully followed, with all trajectory parameters 
well within the limits. 

• Control algorithm design depends only on mass property of RLV providing robustness that 
leads to risk and cost reduction. 

• Reduced development costs is demonstrated by significant reduction of design/tuning 
parameters from 7500 in gain-scheduled PID controller to 36 in TV-SMC. 

• Capability for cost reduction is demonstrated by saving 455 lb of the RCS propellant. 

Test 2. Dispersion case evaluation. 

The goals and implementation conditions of the dispersion case evaluation test are Robustness to 
RLV and environment dispersions are to be demonstrated in 75 dispersion runs using TV- 
SMC/adaptive gain algorithm. 

• The aerodynamic model and actuator model is based on X-33 vehicle data 
The results of dispersion runs are demonstrated for descent only. The bandwidth adaptation in the 
outer loop TV-SMC is shown in Figure 13. The angle of attack, bank angle and sideslip angle 
dispersion evolution is demonstrated in Figure 14. Angle of attack and bank angle tracking errors 
are depicted in Figure 1 5 and demonstrate a reasonable dispersion tracking accuracy. Elevon and 
flap dispersion deflections are shown in figures 16 and 17. They demonstrate acceptable 

performance. The RCS propellant used in descent is shown in Figure 18. The ascent dispersion runs 
are also successful. 
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Summary of the Test2 ResultsAW 75 dispersion trajectories were successfully followed in ascent and 
descent, with all trajectory parameters well within the limits. 

• Capabilities of significant saving of RCS propellant, 340 lb, which amounts to $340,000 
saving if the weight is used for payload, is demonstrated. 

Test 3. Actuator-failed flight evaluation. 

The goals and implementation conditions of the actuator-failed descent flight evaluation test are 

• Demonstrate capability for providing adaptation to non-catastrophic failures and increase 
mission success rate by means of changing bandwidths via TV-SMC. 

• Control commands shall be adaptive to failures and degraded performance. 

• The aerodynamic model and actuator model is based on X-33 vehicle data. 

In Figure 19 inboard elevon deflections are demonstrated. It is clear that the left inboard elevon is in 
a hard on position (actuator failure). Time-varying TV-SMC outer loop gains demonstrate SMC 
gain adaptation in Figure 20. A corresponding angle of attack, bank angle and sideslip angle 
tracking is shown in Figures 21, 22 and 23. The plots demonstrate a high tracking accuracy in 
presence of actuator failure. 


23 












1.00 
0.75 
0.50 
0.25 
0.00 
- 0.25 
- 0.50 
- 0.75 
- 1 .00 
- 1.25 
- 1.50 

200.0 250.0 300.0 350.0 400.0 450.0 

225.0 275.0 325.0 375.0 425.0 475.0 

Fig. 23 Sideslip angle tracking (degrees)/ failure case 
Summary of the Test 3 Results 

• The nominal trajectory was successfully followed in descent with the left inboard elevon in 
hard on position, with all trajectory parameters within the limits. 

• Capability for providing adaptation to non-catastrophic failures and increase mission 
success rate by means of changing bandwidths via TV-SMC/adaptive gain algorithm is 
demonstrated. 

• Control commands are adaptive to failures and degraded performance. 

8. SUMMARY AND CONCLUSIONS 

The time -varying sliding mode controller (TV-SMC) design algorithm is developed for the 2 nd 
generation reusable launch vehicle in ascend and descend modes. In order to maintain stability, the 
bandwidth of the nominal (reduced-order) system is reduced accordingly using a time-varying 
bandwidth PD-eigenstructure assignment technique. The presented TV-SMC fault-tolerant 
technique automatically handles momentary saturations and integrator windup caused by excessive 
disturbances, guidance command or dispersions under normal vehicle conditions. The TV-SMC 
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algorithm has been coded and successfully simulated for the X-33 technology demonstration 
vehicle in ascent and descent modes using high fidelity 6DOF mathematical model. 
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